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ABSTRACT 

Observed protoplan e tary disks consist of a large amount of micrometer-sized particles. 
IDullemond fe Dominikl (|2005D pointed out for the first time the difficulty in explaining the strong 
mid-IR excess of classical T-Tauri stars without any dust-retention mechanisms. Because high rela- 
tive velocities in between micrometer-sized and macroscopic particles exist in protoplanetary disks, we 
present experimental results on the erosion of macroscopic agglomerates consisting of micrometer-sized 
spherical particles via the impact of micrometer-sized particles. We find that after an initial phase, 
in which an impacting particle erodes up to 10 particles of an agglomerate, the impacting particles 
compress the agglomerate's surface, which partly passivates the agglomerates against erosion. Due to 
this effect the erosion halts within our error bars for impact velocities up to ^30 m s -1 . For larger 
velocities, the erosion is reduced by an order of magnitude. This outcome is explained and confirmed 
by a numerical model. In a next step we build an analytical disk model and implement the exper- 
imentally found erosive effect. The model shows that erosion is a strong source of micrometer-sized 
particles in a protoplanetary disk. Finally we use the stationary soluti on of this model to explain the 
amount of micrometer-sized particles in observational infrared data of Fu rlan et al.l (|2006[) . 

Subject headings: accretion, accretion disks — agglomerates — dust — erosion — growth — impact — 
solar system: formation 



1. INTRODUCTION 

Early theoretical investigations concerning planet for- 
mation followed two paradigms. The first paradigm 
considers a laminar nebula in which dust particles co- 
agulate as a result of Brownian motion and differen- 
tial drift velocities. Most particles arriving at the mid- 
plane at 1 AU reach sizes of 1 - 10 cm within 1500 - 
3000 orbits (|Weidenschillindll980HNakagawa et al.lll981l 
1986). This settling leads to a dust-dominated sub 
disk, whi ch keeps a finite thickness o wing to shear tur- 
bulence (iSchrapler fe Hennind 12004 ICuzzi et al.l 119931 : 
lWeidenschillingll2006L 11988111984 11980T ). The shear tur- 
bulence is generated by the difference in orbital veloc- 
ity between the pressure-s upported gas disk and the 
dust-dominated sub-disk (Schrapler & Hcnning 2004; 
ICuzzi et al.lll993T) . In the sub-disk, particles ca n grow to 
large agglomerates (Schrapl er fc Hen ning 2004), whereas 
below and above the sub-disk, a strong outward-directed 
wind is generated. Th is wind, which has typical velocities 
of 30 m s" 1 at 1AU (| Weidenschillind [20061: ICuzzi et al.l 
[1991 carries along small particles that are perfectly cou- 
pled to the gas. 

The second paradigm considers a fully turbulent neb- 
ula. The dust within the disk coagulates due to ran- 
dom velocities provided by the turbu lent gas motion 
(|Brauer et al.ll2008HSchmitt et al.ll9"97b . The strong tur- 
bulent motion in the gas prevents the formation of a dust 
sub-disk. After the particles have coagulated to several 
decimeters in size, they begin to migrate inwards with 
velocities up to 50 m s -1 , owing to friction with th e sub- 
Keplerian gas disk (jWeiden schilling & Cuzzi 1991). 

More recently it was found that the truth may lie in 
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a combination of both models (|Gammid il996: St epinskil 
ll999HCiesiall2007| ). The MHD simulations showed, that 
in a fully turbulent Kelvin-Helmholtz driven disk, tur- 
bulence free "dead zones" appear at places where dust 
prevents the ionization of gas. In this modern paradigm 
of protoplanetary dust growth, fast migration of large 
dust aggregates as well as fast winds (above/below the 
sub-disk) occur at different regions in the nebula. If 
large dust aggregates and micrometer-sized particles co- 
exist in the nebula, then all of the above scenarios lead 
to high-speed collisions between them. Thus, to under- 
stand planetesimal formation under the various scenarios 
described above, it is essential to understand the col- 
lision physics of the dust aggregates embedded in the 
protoplanetary disk. In a series of experim ental papers, 
which are listed in iBlum fc Wurml ([20081 ). the growth 
models of coagulati on by Brownian motion and rel ative 
drift motion s used inlWeidcnschilling fc Cuzzil ([1993D and 
ICuzzi et al.l ([1993D were basica lly confirmed and refined. 
A recent dust-growth model bv lZsom "etTaTl (pOloh . which 
is based upon an elab orate dust-aggregate collision model 
([Guttler et al.ll2010| ). shows that the maximum dust ag- 
gregate size in a minimum-mass solar nebula model is 
around 1 cm, due to the overwhelming effect of bouncing 
in aggregate- aggregate collisions. 

Thus, realistic dust-growth models have two problems: 
(1) They fail to describe the formation of meter-sized 
(and larger) particles, and (2) they fail to explain the 
observed large amount of micrometer-sized particles in 
proto planetary disks without a du st retention mecha- 
nism (jDullemond fc Dominikl I2005D . While the latter 
is a real shortcoming of the models, the former might 
be an indication of the non-existence or short survival 
lifetimes of large protoplanetary dust aggregates. To 
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overcome the d e pletion of smal l grains lOkuzumil (120091 
lOkuzumi et~aT1 (|2010al ) and lOkuzumi et al.l (120101: 



found that in case they are charged their coagulation 
is supressed. 

In this article, we hypothesize that dust aggregates of 
considerable size (e.g. meter size) exist and we follow 
the evolution of these bodies in the protoplanetary disk. 
We will experimentally show that high-velocity impacts 
of micrometer-sized grains into large dust agglomerates 
are a strong source of further micrometer-sized particles, 
which considerably limits the size of the largest surviving 
bodies. 

The experimental method and the laboratory setup 
to perform collisions between micrometer-sized particles 
and dust aggregates are described in Sect. [2] In Sect. 
G2 the experimental results and their analysis are given, 
along with a numerical model that explains and con- 
firms the experimental results. In Sect. |4l the exper- 
imental outcomes are applied to turbulent protoplan- 
etary a-disks with values of the turbulence parameter 
a = 10~ 2 and a = 10~ 4 , which shows the consequences 
of the erosional process on the population of micrometer- 
sized and decimeter-sized particles at different Kepler 
radii. In this Section, we also compare our findings to 
recent infrar e d obs ervations of protoplanetary disks by 
iFurlan et all (|200l . 

2. EXPERIMENTAL APPROACH 
2.1. Mechanical Part 

To reproduce the collision scenario following the 
paradigms described in Sect. [TJ which preferentially 
leads to collisions between the smallest and the largest 
constituents of the dust size distribution, we developed 
a laboratory experiment in which collisions between 
micrometer-sized dust grains (projectiles) and cm-sized 
porous dust aggregates (targets; simulating the surfaces 
of much larger bodies) can be studied. As an analog dust 
material for the projectiles we use sp herical SiC>2 parti- 
cles w ith a diameter of 1.5 /im (see iBlum fc Schraplen 
(2004)). Thus, our projectile particles have sizes typical 
for protoplanetary dust, consist of a silicate- like material, 
and have a morphology suitable for theoretical modeling. 
The targets were produced by random ballistic deposi- 
tion (RBD) out of the same pa rticles and have a volume 
filling factor of 6 = 0.15 (see IBlum fc Schraplerl l200l 
IBlum et "aD l2006f ). RBD agglomerates have the same 
volume filling factor (and, thus, porosity) as the ballis- 
tic particle-cluster agglomerates which should dominate 
the dust-aggregate structures if aggregates preferentially 
grow by collection of small grains. 

To simulate the projectile-target collisions hypothe- 
sized for the protoplanetary disk, a jet of fast projec- 
tiles is directed to centimeter-sized targets in a vacuum 
environment. The mass loss or gain of the targets is de- 
termined by precision- weighing of the targets before and 
after exposure to a well-determined flux of projectiles. 
The experimental apparatus consists of a dust-dispersion 
cogwheel, a velocity filter and a vacuum chamber, which 
contains the target (see Fig. [IJ . The inside of the appa- 
ratus has a residual gas pressure of <4 Pa to prevent the 
particles from losing kinetic energy due to friction with 
the gas. The surfaces of the inside of the experiment are 
covered with vacuum grease on which the particles in- 



evitably stick upon impact. Thus, only particles on bal- 
listic trajectories from the cogwheel through the velocity 
filter can reach the target. To avoid ejected particles from 
falling back to the surface of the target agglomerate, the 
dust jet is directed in a vertically upward direction. 
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pressure inside 
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Fig. 1. — Schematics of the experimental setup, consisting of 
a vacuum chamber, which houses the target, the velocity filter 
(see upper inset photograph), the cogwheel (see lower inset photo- 
graph), the charge-measurement capacitor, and the measurement 
electronics. 



The particle jet is produced with a fast rotating cog- 
wheel (see Fig. [T]), which also deagglomerates the powder 
dust sample into single monomer grains. The cogwheel 
is fed with dust from 20 dust reservoirs, which are ar- 
ranged on a revolver. A piston presses the dust from 
a chosen reservoir against the fast rotating cogwheel. 
The deagglomerated particles possess up to twice the cir- 
cumference velocity of the cogwheel, because they nearly 
elastically rebound from the cogs. The small amount 
of particles, which are not properly deagglomerated, im- 
pact the cogwheel plastically and possess velocities up 
to the circumference speed of the cogwheel. Due to this 
difference in velocity, it is possible to separate the non- 
deagglomerated particles from the monomer grains by 
a velocity filter. The velocity filter consists of an im- 
peller with 60 blades (see Fig. [T]). The inclination of 
the blades is chosen such that their circumference cov- 
ers 20 degrees between top and bottom of the impeller. 
If the impeller rotates at a constant speed, only parti- 
cles within a preset velocity range can pass through the 
device. Particles that are either too fast or too slow 
hit the leading or the trailing blade of the impeller on 
which they stick, because the blades are covered with 
vacuum grease. We measured the velocity distribution 
of the particles leaving the filter by illuminating the par- 
ticles with a pulsed laser beam and observing their tra- 
jectories by a long-distance-microscope equipped with a 
high-speed camera. We found that the velocity of the 
particles leaving the velocity filter possess a full-width- 
at-half-maximum (FWHM) of 17% of the optimum ve- 
locity given by the rotational speed of the filter. The 
velocity filter can be operated from m s" 1 to 90 m s _1 , 
but the cogwheel's deagglomeration efficiency drops dra- 
matically at rotational speeds slower than 7.5 m s _1 . To 
calibrate the mass flow to the target, a greased target is 
placed inside the apparatus instead of the agglomerate 
target so that impacting particles will completely stick 
there. Thus, the mass arriving at the target area can 
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be measured as a function of the dust mass provided to 
the cogwheel. In the second step of the measurement, 
the agglomerate target is weighed before and after the 
impact of a certain amount of particles. The strength of 
erosion (or growth) of the target can be determined by 
the ratio of the mass difference of the agglomerate before 
and after the exposure to the dust flux and the impacted 
mass. In a typical experiment, the target agglomerate 
is placed in the experimental chamber and the chamber 
is then evacuated. Upon reaching the desired vacuum, 
1 to 3 dust cylinders are pressed against the rotating 
cogwheel. After each cylinder, the camber is carefully 
vented to electrically discharge the target agglomerate. 
Then, the agglomerate is removed from the chamber and 
weighed. After weighing, the target is placed back in- 
side the experimental chamber. In typical experiments 
this procedure is repeated 2-10 times, depending on the 
experiment type. 

2.2. Electrical Part 

As known from previous expe riments 

(jPoppe fc Schraplerl [20051 : iPoppe et all I2000bf >. par- 
ticles in the dust jet produced by a cogwheel are 
heavily charged and possess a power-law charge- velocity 
distribution. To assess the influence of the particle 
charge on our erosion measurements, we measured the 
charge of the impactors and the charge transfer to the 
target agglomerate. This was done using the effect of 
electrical influence. The electrical setup is also shown in 
Fig HI The target agglomerate is placed on one insulated 
plate of a measurement capacitor. If an electric charge 
is deposited on the target agglomerate, an identical 
charge is induced on the capacitor plate, which drives 
an electric current to the opposing capacitor plate. This 
current is amplified and logged by a digital oscilloscope, 
which is connected to a computer. The charge on the 
dust agglomerate is determined by integration of the 
electric current over time. To measure the charge of the 
projectile s and compare our measur ement method with 
those of IPoppe fc Schraplerl (|2005l ). a greased target 
was placed on the measurement capacitor instead of the 
target agglomerate. All impacting particles stick on the 
greased target and keep their charge. The number of 
elementary charges on particles, accelerated to velocities 
of 15, 30, 45 and 60 m s _1 , were measured. A typical 
current profile during an experiment is shown in Fig. [2j 

By weighing the greased target before and after the 
impact of the particles, the charge per particle could be 
determined. The average number of elementary charges 
on each monomer particle accelerated to velocities of 15, 
30, 45 and 60 m s" 1 are 44, 198, 318, and 484 e", respec- 
tiv ely, which is in good agreem ent with the measurement s 
by IPoppe fc Schraplerl (|2053) and IPoppe et all (|2000b| ). 
As a result of the charged dust target, the impactors 
lose kinetic energy before an impact. The calculation of 
the energy loss due to the electric potential was approx- 
imated by assuming the target agglomerate as a charged 
sphere with the radius of the target diameter and the 
charged particles as test charges arriving from infinity. 
These calculations were confirmed by a numerical simu- 
lation of the incoming particles using the exact geome- 
try of the experimental setup with the software package 
COMSOL. In all measured parameters, the velocity loss 
of the impinging projectiles was smaller than 2%. The 
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Fig. 2. — A typical electric current profile when a greased target 
is exposed to a dust jet with 15 m s — 1 impact velocity. The dust 
flow starts at t m 50 s and ends at t ss 300 s. The dust jet is 
generated by pressing one dust cylinder against the cogwheel. 

collision dynamics of the erosion process is not changed, 
because the required pull-off force of the particles is four 
orders of magnitude larger than the electric force close to 
the target surface. Charging of the particles and targets 
is therefore negligible for the following measurements. 

Fig. [3] shows an example of the influence of the projec- 
tile and target charging on the impact energy and veloc- 
ity. In the top left graph of Fig. [3l the electric current 
to the target agglomerate caused by the impinging pro- 
jectiles with a mass of 3.5 x 10~ 15 kg is shown. The 
current drops to zero after a certain number of particles 
have impacted. This is due to the fact that particles are 
deposited to the surface and charges that are separated 
as a cause of the impact are also sitting on the surface. 
Each impact erodes particles from the agglomerate sur- 
face (see below) so that after a certain number of impacts, 
an equilibrium between deposited and separated charges 
(sitting on eroded particles) is reached. In the top right 
graph of Fig. [3]the total charge of the target agglomerate 
is given as a function of the accumulated mass flux of the 
impinging projectiles. Note that this charge is limited by 
the equilibrium discussed before. In the lower left graph 
of Fig. the total energy loss caused by Coulomb re- 
pulsion of an impinging particle is shown. In the lower 
right graph of Fig. [31 the velocity loss of an impinging 
projectile is given, which is negligible compared to the 
initial velocity of 60 m s . 

3. EXPERIMENTAL RESULTS 

3.1. Erosion and Passivation 

We performed the impact experiments described in 
Sect. [2] at impact velocities of 15, 30, 44, and 59 m s , 
respectively. A surprising result is that the erosion effi- 
ciency (i.e. the relative mass loss of the target) depends 
on the impact history of the (initially highly porous) tar- 
get agglomerate, as can be seen in Fig. @]for impacts with 
59 m s^ 1 projectile velocity. The data points show the 
mass loss of the target agglomerate, Am, per imping- 
ing monomer mass, m; = 3.5 x 10 -15 kg, as a function 
of the total mass exposure per unit cross section onto 
the target, ijl. One can see from the saturation of the 
curve that the target agglomerate is passivated against 

erosion for mass exposures fi ~ 2 mg cm~ 2 . This sat- 
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Fig. 3. — Examples of the influence of electric charges on the 
collision behavior of micrometer-sized dust projectiles (with a mass 
of 3.5 X 10 — 15 kg) impacting a porous dust-aggregate target at 60 
m s _1 velocity. Top left: the profile of the electric current over 
the total amount of dust exposure. Top right: the total charge 
accumulated by the target. Lower left: the calculated energy loss 
of a charged projectile carrying 1000 elementary charges. Lower 
right: the calculated velocity loss of a charged projectile carrying 
1000 elementary charges with an initial velocity of 60 m s" . 

uration mass exposure is about 20 times the number of 
particles on the surface layer of an unprocessed agglom- 
erate. Experiments with impact velocities of 15, 30, 44, 
and 59 m s -1 show a similar saturation behavior. The 
error bars in Fig. U denote the measurement uncertain- 
ties of the weighing process. The erosion efficiency at the 
initial and final stage for all measured impact velocities 
is shown in Fig. [5] The velocity error bars denote the 
velocity dispersion due to the width of the velocity filter. 
The error bars in the vertical direction are the standard 
deviation of several measurements with different target 
agglomerates. Fig. [5] shows that the initial erosion ef- 
ficiency scales with impact velocity but is reduced to a 
saturation level of zero for impact velocities lower than 
26 m s _1 . For higher impact velocities, the erosion effi- 
ciency is reduced to about 15% of its initial value. The 
reduction of the erosive effect of the particle impacts goes 
along with a compaction of the surface. Figure [5] shows 
scanning electron microscopy (SEM) images of the tar- 
get agglomerate's surface before and after an exposure 
to dust fluxes at 15, 30, 44, and 59 m s _1 , respectively. 
An obvious change in surface structure occurs between 
30 and 44 m s" 1 . 

3.2. Theoretical Considerations 

3.2.1. Analytical Part 

A key to understanding the erosion behavior describe d 
above is given in the work of IKonstandopouiosI ((2000). 
He numerically simulated the formation of agglomerates 
by the RBD process and discovered that the coordina- 
tion number n (i.e. the number of next neighbors) in an 
agglomerate increases with increasing impact velocity of 
the depositing particles, until a maximum coordination 
number of n max = 6 is reached. Note that the max- 
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Impact History (mg/cm 2 ) 

Fig. 4. — Erosion efficiency of monomer projectiles impinging a 
target agglomerate at 59 m s — 1 . The saturation of the erosion 

efficiency for fj, ~ 2 mg cm -2 is clearly visible. The error bars 
result from the measurement uncertainty of the weighing process. 
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Fig. 5. — Erosion efficiency of target agglomerates for impact ve- 
locities of 15, 30, 44, and 59 m s" 1 , respectively. The upper data 
points with error bars were recorded at the onset of the experiment, 
in which the target agglomerate was very porous and not yet passi- 
vated, whereas the lower data points with error bars show a reduced 
erosion efficiency after the target-agglomerate's surface was almost 
passivated against erosion. The velocity error bars denote the ve- 
locity distribution stemming from the velocity filter. The vertical 
error bars denote one standard deviation of several measurements 
with different target agglomerates. The solid curve follows Eq. 151 
and i s the result of a numerical impact model described in Sect. 
l3~T2l 

imum coordination number found by IKo nstando pouiosI 
(2000) has a gaussian d istribution with a FWHM of 2. 
IKonstandopouiosI ([2000D found that a relation fits his nu- 
merical results which reduces the collision of a single par- 
ticle with a particle bound in a large agglomerate to a 
collision of a single particle with another single particle 
of an effective mass given by 

m e g — m(l + Cn), (1) 

where C is a dimensionless deposit "rigidity" parameter 
in the range (l,oo), m is the mass of the particle bound 
to the agglomerate and n its local coordination number. 
Using Eq. [T] and applying momentum and energy equa- 
tions we calculated the energy transferred to an effective 
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Before Impact 



After Impact of 15 m/s Monomers 




After Impact of 30 m/s Monomers After Impact of 44 m/s Monomers 




Fig. 6. — SEM images of the surface of a target agglomerates 
before and after the passivation by dust-particle impacts at veloc- 
ities of 15, 30, 44, and 59 m s , respectively. The scale bars at 
the bottom left of each SEM image are 20 fim for the bottom-right 
image - and 10 ^m for the other images. 



particle to 



1 + Cn 
'(2 + Cn)< 



(2) 



with E imp denoting the kinetic energy of the impactor. 
The energy to release a particle from an agglomerate is 
given by 



E r = E c n, 



(3) 



where E c is the break-up energy of a single particle- 
particle contact. An empirical relation for the capture 
threshold velocity v s t can be deduced from the experi- 
ments bv lPoppe et al.l ()2000af ). who showed that 



v, t = 0.92 



/im 



-0.52 



m s 



(4) 



with ao being the monomer radius. Based on this re- 
lation, we can derive the break-up energy of a particle 
contact to be 



E = 1 5 



m v 2 st , 



(5) 



(Blum k Wurmll200(fl with m being the mass of the im- 
pacting monomer. It is obvious from Fig. |5]that the tar- 
get agglomerate gets compacted by the impinging dust 
projectiles. An increased compression, i.e. an increased 



volume fillin g factor tp results in a n increased coordina- 
tion number (Lagcmaat ct al. 200lJ) . The increased coor- 
dination number increases the release energy of a particle 
and is a mechanism that passivates the agglomerate. 

3.2.2. Numerical Description of an Unpassivated Target 
Agglomerate 

iDominik fc Tielensl (fl997l ) and lWada etaLl pOOl pre- 
dict an erosion effect proportional to the impact energy, 
thus proportional to the squared impact velocity, which is 
in contradiction to our experimental results (see Fig. [5J. 
To validate our measurements by theory, we developed 
a numerical model that applies the collision model of 
iChokshi et al.l (119931) wit h the adaptation recommended 
by iBlum fc Wurml (|2000ft . In our model, the energy loss 
of a repelling particle equals its binding energy and t he 
effective mass term (Eq. [1]) of lKonstandopoulosI (|2000f ) is 
used. The effective mass term r equires a rigidity param - 
eter C. For fluffy agglomerates iKonstandopoulosj (2000) 
chose it close to 1 to get realistic results. Therefore we 
also chose a value 1 for the rigidity parameter. In our 
model, all particles of an agglomerate have an effective 
mass corresponding to their coordination number and 
thus only part of the impactor 's kinetic energy is trans- 
ferred to the respective agglomerate monomer. The re- 
maining kinetic energy of the impactor after its impact 
is either used in a further impact or is lost if the particle 
escapes from the agglomerate. The escape probability 
S of the impactor and the ejected particles is calculated 
from the surface texture of a num erically grown RBD 
agglomerate (see IBlum et al.ll2006[ ). This calculation is 
done in three steps. In the first step, all particles that 
can be reached by impinging particles from a direction 
perpendicular to the (mean) target surface are identified 
and defined as surface particles (see Fig. [7]). Depending 
on the fraction of the particles on the aggregate's surface 
that can be reached by impinging particles of the same 
size, a corresponding accessability is calculated. In a sec- 
ond step, the angles of the surface particles are calculated 
under which they can escape without impacting again 
with the target agglomerate after they are released by 
the original projectile. The same is done for the reflected 
impactor. For simplicity this procedure is not applied to 
all particles and neighbor particles on the agglomerate 
surface. The problem was reduced to a 2D-problem. 

Only particles on a given line on the surface are taken 
into account, and to calculate the required angles, only 
particles are used that are found in one direction of that 
line. It is therefore assumed that the target surface is 
isotropic. The smallest of the angles from a given particle 
on the line to its neighbors in one direction is the maxi- 
mum escape angle (see Fig. [7]). By repeating this proce- 
dure for a large number of particles on the defined line 
on the agglomerate's surface and by averaging over the 
occurring angles, the mean escape angle is determined. 

The probability to impact a further particle is calcu- 
lated by assuming that this is the case if the reflection 
angle is larger than the mean escape angle. The prob- 
ability of a particle to fulfill this condtion is found by 
taking into account that all particles are spherical. The 
reflection angle of a particle impinging perpendicular to 
the agglomerates surface is connected with the impact 
distance from the target particle's axis ri in the agglom- 
erate (see Fig. [7]). 
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Fig. 7. — Left picture: An impacting particle (white) is shown 
that is reflected off an agglomerates surface at an angle tp. The 
black particles are surface particles that can be reached by the 
impinging particle. The grey particles are bulk particles of the 
agglomerate that cannot be reached. The angle ip corresponds 
to the probability that the impactor can impact a second time. 
Right picture: Relation between the maximum escape angle and 
the probability of a second impact of the impinging particle r2r _1 . 

The correspondence between the reflection angle and 
the probability of an reflection angle larger than <j> is 
given by 



'"2 



which corresponds to 



2 = 1 — sin — 
2 



(6) 



(7) 



(see Fig. [7]). These probabilities are multiplied with 
the corresponding accessibility of that given particle and 
averaged over all considered particles. We found that for 
RBD agglomerates 3 = 0.39. 

If the energy transferred to the agglomerate target ex- 
ceeds the binding energy of the monomer grain, it will 
be ejected and will impact another agglomerate monomer 
with a probability of 1 — S. This corresponds to the fact 
that the mean coordination number of the (initial) target 
agglomerate is 2, which means that it is an agglomerate 
of particle chains. It is very unlikely for an impactor to 
hit the end of a chain in a direction parallel to the chain. 
In all the other impacts the hit particle can be directly 
ejected. This is not the case for a compacted agglomerate 
with a larger mean coordination number. 

The contact energy E r of a monomer p articles re- 
quir ed for the modified [Chokshi et al.1 (|1993j) model (see 
alsolDominik fc Tielenslll997[ ) is taken from lPoppe et all 
^OOOajj - The coordination number distribution of the 
particles at the agglomerate's surface and the surface 
texture of the agglomerate are also c alculated from th e 
numerically grown RBD agglomerate (|Blum 
Our numerical erosion model is realized as a statisti- 
cal code with recursive elements. This means that all 
possible impacts (concerning the effective mass of the 
target particles) of an impinging particle are calculated, 
multiplied with the probability to release a particle and 
summed up. All possible impacts of the reflected im- 



pactor and the ejected particles are multiplied with the 
probability of a further impact and the probability of 
the occurrence of the effective masses of the targets and 
are also summed if their impact energy is large enough 
to release a particle. This is done in a recursive way 
until the particle energy of a given branch of the simu- 
lation is lower than E c which then leads to sticking. In 
other words, to get the total number of particles that 
are ejected by the impact of a given particle, the proba- 
bilities of all possibilities to eject a particle are summed 
up, where the possibility to eject a certain particle is 
the product of all possibilities of the preceding path of 
the impactor. The outcome of the impact simulation on 
the unprocessed agglomerate is shown as a solid curve in 
Fig. [5l It can be fitted for impact velocities larger than 
15 m s _1 by a squareroot function 



Am 



1.25 v 



imp 



1 m s 



7.5. 



(8) 



In the region lower than 15 m s _1 sticking can occur, 
which considerably reduces the erosion efficiency. For 
velocities below 2.7 m s _1 , sticking dominates the impact 
process, which leads to a net growth of the agglomerate 
i.e. Am/rrii < 0. At velocities lower than 1 m s _1 all 
impactors stick i.e. Am/mi = — 1. 

3.2.3. Numerical Description of a Passivated Target 
Agglomerate 

In a passivated agglomerate, the projectiles impact 
perpendicular to the surface of interconnected particles. 
The ejection of such a particle requires more energy 
than for the loosely-bound monomers in the RBD tar- 
gets. As no numerical agglomerate for passivated tar- 
gets is available, we have to rely on the structural infor- 
mation provided by the SEM pictures, which are two- 
dimensional. Therefore it is not possible to exactly mea- 
sure the coordination number of a passivated agglomer- 
ate. We also choose C=l for t his simulation. Th e coor - 
dination number of 6 found bv iKonstandopou los (2000) 
leads to a too perfect passivation at nearly all veloci- 
ties. Koji Wada (personal communication 2009) found in 
molecular-dynamics calculations that agglomerates that 
are processed by collisions tend toward a mean coordina- 
tion number of 4. Therefore we chose a mean coordina- 
tion number distribution on the agglomerate surface of 3 
with a gaussian distribution of a FWHM of 2. The passi- 
vated agglomerates show a steep valley-hill-structure (see 
Fig. [S] bottom right). 

SEM images of the type shown in Fig. [5] were used to 
estimate the probability S p that an ejected or reflected 
particle will again impact at the agglomerate. In a passi- 
vated agglomerate, the monomers bound in the agglom- 
erate cannot be directly ejected as it was the case with 
the unpassivated agglomerates, which mainly consists of 
particle chains. After an impactor hits a monomer bound 
in a passivated agglomerate, the bound monomer itself 
impacts its up to 5 neighbors. If its kinetic energy after 
this multiple impact is higher than the binding energy 
to its neighbors, the monomer is ejected. In our nu- 
merical code, all possible impacts with respect to the 
coordination number distribution were calculated and 
multiplied with the estimated occurrence of each coor- 
dination number. Any target monomer that is hit by 
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any other monomer particle (projectile or ejected target 
monomer) can again impact with its neighbor particles in 
a recursive way. All ejected particles can again impact on 
the agglomerate with the probability of H p . Please note 
that because of the lacking information on the internal 
structure of the passivated agglomerate, this numerical 
method can only be a rough estimate. As a result of our 
crude model for passivated targets (see Fig. [8]), it was 
found that the order of magnitude of the agglomerate 
passivation against erosion can be reproduced by the nu- 
merical model. However, due to the crude assumptions of 
our model, it is not able to reproduce the detailed charac- 
teristics of the erosion efficiency as a function of impact 
velocity. Thus, for the following application (Sect. [4|, 
the measured values are used via a linear interpolation 
function of the measured data 



Am 



u imp 



(9) 



with a — 0.0155 s m 1 and b — —0.4 in the velocity range 
26m s" 1 < v < 60m s _1 . 
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Fig. 8. — Erosion efficiency as a function of impact velocity for 
a passivated agglomerate target. The crosses denote the experi- 
mental values with error bars, the solid line shows the result of our 
numerical model. The dot dot dashed lines mark the zero erosion 
line and the fit for the velocity dependent erosion given in Eq. 
and used in Sect. [4] 



4. APPLICATION TO A PROTOPLANETARY DISK 

4.1. Grain Population in a Solar Nebula 

To estimate the influence of the erosion process on the 
dust-aggregate distribution in a protoplanetary nebula, 
a simple analytical model based on the timescales of the 
erosion and coagulation processes will be developed and 
compared to observational findings in the following sec- 
tion. For this we use a standard minimum mass proto- 
planetary nebula with a disk mas s of 0.0425 M^. The 
analyt ical disk model is taken from[Schrapler fc Pfenning! 
(2004) (their Appendix Al). This leads to a surface mass 
density of the at 1 AU of 1.7 x 10 4 kg m -2 and a gas 
density in the midplane at 1AU of 1.4 x 10~ 6 kg m~ 3 
and 280 K at 1 AU. The dust mass is 1% of the gas 
mass. The relative v e locitie s of the particles are taken 
from [Wcidcnschillingl (|1977t ) (his Figures 3 and 5) and 
approximated and extrapolated in Eqs. 1111 1551 and I56|) . 
The number of large particles n m is calculated for a given 



Kepler radius by assuming that 70% of the total mass 
is in the large dust agglomerates with a size given in 
Figs. [10] and [TTJ (bottom graph, solid line). The total 
mass ratio of 70% for the large du st aggregates is taken 
from INomura fe Nakagawal ([20061 ) (their Fig. 5) after 
one coagulation timescale, the radii of the large dust 
aggregates are taken from INomura fc Nakagawal (2006) 
(their Fig. 6, dash-dotted line). A increase of the total 
mass in the large agglomerates to 90% of the total mass, 
which would be the ca s e afte r 2 coagulation timescales in 
INomura fc Nakagawal (2006) (their Fig. 5) leads to mi- 
nor changes of our results that wuold be hardly visible 
in Figs. [TU] and [TTJ We ignored icy grains at outer disk 
radii. The threshold impact velocity for erosion should 
be considerabely higer for icy particles. However the pre- 
dominent part of the erosion takes place at Kepler radii 
below 2AU. 

4.1.1. The Erosive Formation of Monomer Dust Grains and 
Their Subsequent Coagulation 

Most of the coa gulation models (e.g. 
INomura fc Nakagawal 120061) show that micrometer- 
sized particles vanish from the nebula on a very short 
timescale, due to their effective coagulation. Therefore, 
it is important to calculate whether the source of small 
dust grains released by the erosion of large dust aggre- 
gates, as described in the previous sections, is strong 
enough to compensate the loss of the dust particles by 
coagulation. As we have shown above, a small initial 
number of micrometer-sized dust grains is sufficient to 
start the erosion process of aggregates with sizes of 
decimeter and above. The particle production rate in a 
differential volume element dU = 2nRHdr (see Fig. [§]) 
is given by 



d(dAT) 
dt 



Pa n v imp f dA 



(10) 



where the relative velocity between the micrometer- 
sized particles a n d the large agglomerates is taken from 
Weidcnschilling] ([1977D and approximated for Kepler 
radii larger than 0.2AU to 



50— RfAU] 

s 



(11) 



Here, f(vi mp ) is the number of monomer grains ejected 
by an impacting particle, po is the number density of 
micrometer-sized particles in the midplane of a proto- 
planetary disk, and n is the number of large dust aggre- 
gates at a disk area Ary — 2irRr, where it is implicitly 
assumed that most of these large aggregates are located 
close to the midplane. This area element is chosen be- 
cause all particles located in there cut the differential 
volume element dU with their front side in the velocity 
direction. Therefore the number n can be expressed by 
the surface density 5 of large agglomerates 



n = 52irRr. 



(12) 



The mean differential area with which the particles in- 
side Ab cut dV, is given by 



dA r , 



dr 

:27T — 

r 



s ds. 



(13) 
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Fig. 9. — Upper picture: The differential volume element dV, 
which is a torus in the protoplanetary nebula at the orbital distance 
R with a height H and a thickness of dr. Lower picture: A large 
particle with a radius r at a radial distance s from dV in the 
midplane of the disk. As s < r an area of intersection forms as a 
ring with a radius s and a thickness dr. 



With 



we get 



h 2 and ds 



hdh 



\/r 2 — h 2 



dA r , 



rdr, 



(14) 



(15) 



assuming that all particles are homogeneously dis- 
tributed in 4o. Here, h is the distance of the center 
of a given particle to dV (see Fig. [9]). 

The number of micrometer-sized particles produced 
within a given radius element 2irRdr is diV, which will in 
the following be expressed using the corresponding par- 
ticle number density in the midplane po. The differential 
surface density of the particle number diV at a given ra- 
dius element 2nRdr is 



_dN_ 

2irRdr ' 



The integration 



E = 



p(z)dz 



of the vertical number density structure 

^2 \ 



P{z) 



p exp 



2H 2 ) 

in the protoplanetary nebula results in 

E 

Po 



HV2n 
Inserting Eq. [16] yields 

dN = p a (2Tr)?RHdr 



and 



df dt v ' 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



Inserting Eqs. HUrTJl and [15] into Eq. [TU] results in 



dpo 
dt 



Po- 



- 'imp 



V2H 



(22) 



which is the differential equation for the dust number 
density increase of micrometer-sized particles without co- 
agulation. 

In the following step we will introduce the coagula- 
tion of the pm-sized particles. For simplicity and as a 
lower limit to the erosion efficiency we assume that all 
micrometer-sized particles that coagulate to dimers or 
larger aggregates are lost to the erosion process. This 
assumption leads to the equation 



d P 2 
- = -ovsp 



(23) 



( Blum 2004), where we assumed that the main source for 
the inter-particle collisions is the Brownian motion with 
a mean collision velocity vb and a is the cross section 
in the collision of two micrometer-sized particles and p is 
given by Eq. [18] Because /j,m-sized particles are perfectly 
coupled to the disk gas, turbulence effects are negligible. 
To calculate the coagulation particle losses we have to 
take into account that the coagulation efficiency depends 
on p and therefore on the disk scale height H . Because 
dust particles of /im size are mixed over H on a timescale 
of an orbital period it is possible to integr ate both sides 
of Eg . [23l over the disk height following iBrauer et al.l 
( 2008) (in their Appendix B) which results in 



dE 
~dl 



-\firHcrvBpo- 



(24) 



With Eq. [19] we obtain an equation for the midplane 
number density of monomers 



dpo _ mbpI 
dt V2 



(25) 



Performing the coagulation feedback by adding Eq.l2~5lto 
the righthand side of Eq. [2H results in 

dpo 
dt 

with 

r 2 Sv imp f^ 



: Qp - ®pl 



(26) 



V2H 



<JVB 
V2 



whose solution is 



Po(t) =poo 



6 exp(et) 



e + p o*(exp9<- 1) 



(27) 
(28) 

(29) 



Here, poo is the initial number density of micrometer- 
sized dust. The half- value time of the above equation, 



( Aio$ 



e 



0.59 



(30) 



seems to be dependent on poo, but for a typical pro- 
toplanetary nebula ^> poo^j which is valid for any 
reasonable poo- Therefore, the half- value time is 



t 2 



e" ln2 - 



(31) 
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For a minimum mass solar nebula at 1 AU t% is about 
2 years. Because H is roughly proportional to R, fa 
approximately scales with the Kepler radius squared. As 
we will see later, it is of the same order as the timescale 
for erosional radius loss of large particles. Therefore, the 
stationary solution of Eq. [55] is used in the following as 
a zero-order approximation, i.e. 



POs 



u tmp 



Havr 



(32) 



4.1.2. Timescale of the Erosive Radius Loss of Large Dust 

Aggregates 

The number of the monomers eroded from a large dust 
aggregate per unit time is given by 



df 



= — r 7Tiv 



ofPOs 



(33) 



The relation between the radius of a large particle r and 
the number of monomer grains e inside the aggregate is 



(34) 



Here, I is the filling factor of the large dust aggregate and 
Vfj, denotes the volume of a micrometer-sized particle. 
Differentiation by t gives 

de At 2 ixl dr 
dt~ V,, dt' 



(35) 



Combining Eqs. 1331 and 1331 results in 



dr 
df 



"77 : " ~^ v impfPQsV fl . 



(36) 



The value 



imp 



do es not depend on r for particle sizes 



larger than 0.3 m in I Weidenschillind (|1977t) (his Fig. 2). 
Due to an expected filling factor of 0.4 for macroscopic 
protoplanetary dust aggregates and a bulk density of 
2000 kg m~ 3 (see (|Zsom et alJ 120101 ) ). which leads to 
an aggregate density of 800 kg m~ 3 , w e re-ca lculated 
the relative velocities of iWeidenschillingl (|1977f ) for ag- 
gregates of 40% filling factor. For these dust aggregates, 
the velocity does not depend on the aggregate size for 
radii larger than 0.6 m at 1AU, which is in zeroth order 
approximation also the case for larger Kepler radii, with 
ro=lm. Inserting Eq. [32] into Eg. 1361 gives 



dr 
dt 

whose solution is 




AIHcjvb 



AIHcfvb 



The corresponding half value time of r is 

AlHavB 



ti = 



(37) 



(38) 



(39) 



The erosion timescale as a function of the Kepler radius 
is given in Figs. [TU] and [11] (top graphs, dashed-dotted 
lines) for turbulence parameters a = 10 -2 and a = 10~ 4 , 
respectively. The values for / for different velocities were 



taken from the data points of Fig. [5] where it was as- 
sumed that the curves in between the data points are 
straight lines. 
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Fig. 10. — Protoplanetary disk with a = 10 — 2 . Top graph: co- 
agulation timescale (solid line), erosional timescale (dashed-dotted 
line), radial turbulent transport (dotted line). Second graph: num- 
ber density ol micrometer-sized particles at erosion timescale larger 
than the coagulation timescale (solid line), number density of 
micrometer-sized particles at erosion timescale smaller than the 
coagulation timescale (dashed-dotted line). Third graph: ratio 
of the total mass of micrometer-sized particles to the total dust 
mass at a given Kepler radius. Bottom graph: Particle size in 
which 70% of the total solid disk mass is contained (also maxi- 
mum particle size) for er osion timescales larger than t he coagula- 
tion timescale taken from Nomura &: Nakagawa (2006) (triangles) 
and their power-law interpolation (solid line). Maximum particle 
size for erosion timescales smaller than the coagulation timescale 
(dashed-dotted line). 



4.1.3. Passivation Timescale 

To decide whether the erosion efficiencies of the unpas- 
sivated or the passivated agglomerates should be used to 
describe the erosion processes in the solar nebula (see 
Fig. [5]), the passivation timescale needs to be estimated. 
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Fig. 11. — Protoplanetary disk with a = 10 — 4 . Top graph: co- 
agulation timescale (solid line), erosional timescale (dashed-dottcd 
line), radial turbulent transport (dotted line). Second graph: num- 
ber density of micrometer-sized particles at erosion timescale larger 
than the coagulation timescale (solid line), number density of 
micrometer-sized particles at erosion timescale smaller than the 
coagulation timescale (dashed-dotted line). Third graph: ratio 
of the total mass of micrometer-sized particles to the total dust 
mass at a given Kepler radius. Bottom graph: Particle size in 
which 70% of the total solid disk mass is contained (also maxi- 
mum particle size) for er osion timescales larger than t he coagula- 
tion timescale taken from Nomura & Nakagawa (2006) (triangles) 
and their power-law interpolation (solid line). Maximum particle 
size for erosion timescales smaller than the coagulation timescale 
(dashed-dotted line). 

The number of particles on the surface of an agglomerate 
is given by 



N. 



surf ~ 2-nr 2 



(40) 



where r pm is the radius of a monomer. The rate at which 
monomer particles impact on a dust aggregate's surface 
is given by 



dN rl 



dt 



■-por 2 TTv imp f. 



(41) 



We found in our experiments (e.g. Fig. 2]) that an ag- 
glomerate is passivated after about 10 times the number 
of surface monomers have impacted on the agglomerated. 
Thus, the passivation time is 



, dN coU 



tpass — ^ 



(42) 



This passivation timescale t pass is smaller than a day 
for all reasonable parameters and orbital radii in the 
disk. Therefore, passivated agglomerates are dominat- 
ing the erosion process so that we use the experimental 
erosion efficiency from the passivated dust samples (dot- 
dot-dashed line in Fig. [5] and Eq. 

4.1.4. Coagulation Timescale 

For the coagulation timescales we u se the results ob- 
tained by INomura fc Nakagawa! (|2006f ) . Their model is 
based on the assumption that compact particles (with a 
filling factor of 100%) that grow to a certain size are 
removed from the coa gulation pro c ess by radial drift . 
In other models (e.g. IB rural 120041: iBrauer et all 120081: 
iZsom et all 120101) , the coagulation is counteracted by 
fragmentation and bouncing of colliding particles, which 
is not considered here. As macroscopic protoplanetary 
dust aggregates possess fil ling factors of about 4 %, we 
will adapt the model by INomura k, Nakagawal ()2006l ) 
to our non- compact dust aggrega t es. T he coagulation 
timescale of INomura fc Nakagawal (|2006l ) is given by 



Tk 



(43) 



with p m , A m , and v r being the number density of dust 
aggregates with mass m, their collision cross sections, 
and the relative velocities, respectively. 

The coagulation timescale Tk depends on the filling fac- 
tor I via the cross section 



A m ll 



which leads to 



Tk 



J3. 



1 



(44) 



(45) 



where A m n is the cross section of a particle in case I = 1 
and p m is the number density of these particles. The 
relative velocity is approximated for simplicity by the 
velocity caused by turbulence, v r = v tur b, and is assumed 
to be independent of the pa r ticle s ize and filling factor. 

In INomura fc Nakagawal (|2006f) coagulation is bal- 
anced by the radial drift of the particles with a velocity 



Vd = 9dTf, 



(46) 



where gd is the local acceleration caused by the sub- 
Keplerian motion of the particle due to gas drag. The 
friction time of a given particle is 



9 p g Vr, 



(47) 



with p p , p g , rd, and v m being the mass density of the 
dust particle (proportional to I), the mass density of 
the gas, the maximum radius particles can grow, and 
the mean molecular v elocity of the gas, respectively (see 
iWeidenschillingl Il977l ). Here, p p depends on the bulk 
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density pb of the dust monomers and the filling factor of 
the dust aggregate through 



P P = Ipb- 
The drift timescale is given by 



Td 



Sd 

Vd ' 



(48) 



(49) 



where Sd is an arbitrary drift length. As the only depen- 
dency of Eqs. and 2H]on the volume filling factor 
I is in p p , Td is thus proportional to I . This means 
that the coagulation timescale scales with Z§, whereas 

the drift timescale scale s with I -1 . 

To use the results of iNomura fc Nakagawa (2006|), a 
new maximum particle radius (dependent on I) has 
to be found with which coagulation and radial drift still 
balance each other. The ratio of the two timescales, 



2l° 3 gdPbrl 



Tk(l) 

Td(l) 9pmA m iVturbPgVmSd' 



(50) 



is proportional to 1^ and proportional to r 2 .. This means 
the coagulation with a filling factor of I is still bal- 
anced by the radial drift if is multiplied by a fac- 
tor l~~s . Therefore, the maximum particle radius at 
1AU for a volume filling factor of I = 0.40 is about 0.5 
m. Fr om these considerations, it follows that the re- 
sults of lNomura fe N akagawa ( 20061) are still valid for our 
lower filling factors if the maximum size of their parti- 
cles is increased by a factor l~e . The original coagulation 
timescale will then be changed by l~3 . However, the rel- 
ative size distribution derived by INomura fc Nakagawa 
(2006) will remain uncha nged. The original coagula - 
tion timescales are found in INomura fc Nakagawa! (2006) 
(their Fig. 5), which is the time at which 70% of the total 
dust mass is in particles with a critical size at which they 
are lost from the coagulation process by fast radial drift. 
The dependency of the coagulation timescale on the disk 
radius is — 70 years x 1$ for a = 10~ 2 and = 500 
years x It for a = 10~ 4 , respectively. We interpolated 
the timescales for the other radii assuming a power law. 
The resulting coagulation timescale as a function of the 
Kepler radius is shown in Figs. [10] and [TTJ as solid lines 
in the top graph. 

4.1.5. Timescale of Radial Turbulent Transport of 
Micrometer-Sized Dust 

The timescale of turbulent transport T tur b of 
micrometer-sized dust is 



7~turb 



R 2 

Vturb 



(51) 



The turbulent gas viscosity is given by 

vturb = ac s H, (52) 
IShakura fc: Sunvae v (1973), with the sound speed 



hk B T 



(53) 



where 7 is the adiabatic exponent of H2, fcs is Boltz- 
mann's constant, and m g is the mass of an H2 gas 



molecule. The temperature is assumed to be dependent 
on the Kepler Radius Rk via 



T = 280 K 



\1AU 



(54) 



(see lSchrapler fc Henning||2004D . The timescale of turbu- 
lent particle transport as a function of the Kepler radius 
for a = 10 -2 and a = 10~ 4 is shown in Figs. IT01 andim 
respectively (top graphs, dotted line). It is much longer 
than the coagulation and the erosion timescale. There- 
fore the effect of turbulent particle transport is negligible 
for our considerations. 

4.1.6. Particle Size Limit Caused by Radial Drift 

The model of lNomura fc N akagawa] (|2006l ) is based on 
the assumption that compact particles that grow to a 
certain size a re removed from the coagulat ion process by 
radial drift. INomura fc Nakagawa! (|2006l ) find in their 
Eq. 24 and the sentence below that these particles can 
grow close to the midplane to 32 cm, 80 cm, and 6 cm at 
1 AU, 5.2 AU, and 30 AU, respectively for a turbulence 
strength of a = 10~ 2 and a = 10~ 4 . The assumption of 
compact particles should now be dropped in favor of a re- 
alistic assessment of the porosity of the dust aggregates. 
We have seen above that fluffy dust aggregates get com- 
pacted by the continuous bombardment of micrometer- 
sized dust particles. Our model for the erosion efficiency 
(see Fig. [8]), which assumes a mean coordination number 
of 4, fits reasonably well. Such a coordination number is 
rel ated to a volum e fillin g factor by the relation shown 
in lLagemaat et ail (|2001[) . From their Fig. 4, we can 
conclude that the corresponding (surface) filling factor 
of the aggregates was increased to ( « 0.4. To a simi- 
lar result fo r the ( bulk) filling factor, we get if we follow 
iZsom et al.1 ()2010[ ) who calculate that macroscopic dust 
aggregates should have volume filling factors I ss 0.4, due 
to low-velocity bouncing collisions. Thus, we assume in 
the following a volume filling factor of I = 0.4 for large 
dust aggregates in PP Ds. For this, we tran s forme d the 
dust-aggregate sizes of lNomura fc Nakagawal (|2006h with 
the factor calculated in Sect. 14.1.41 These trans- 
formed sizes are plotted as triangles in Figs. [TO] and [TTJ 
(lower graphs). In between the data points, we assumed 
power laws, which are shown as solid lines in the lower 
graphs of Figs. [TOI and [TO 

4.1.7. Particle Size Limit Caused by Erosion 

In INomura fc Nakagawal (|2006l ). the coagulation of 
particles is limited due to the fact that large particles are 
lost by the drift toward the central star. However, in our 
case and at Kepler radii R < 0.2 AU (a = 10 -2 ) and R < 
0.3 AU (a = 10~ 4 ), the erosion timescale is much shorter 
than the coagulation timescale (see Figs. ITOl and ITT1 top 
graphs). Therefore, the particles can only grow until the 
onset of erosion at drift velocities of 26 m s _1 (see Fig. 
[8]) . The dependency of the impac t velocity, Vj mr >i on the 
particle radius is estimated using iWeidenschillingl (]1977fl 
(his Fig. 3) for a particle density of 1000kg m~ 3 . The 
impact velocity can be approximated by 



26 



s V 0.21m 



(55) 
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The dependency of the impac t velocity, Vimn, on the 
Kepler radius is estimated using IWridcnschillmg] ([1977D 
(his Fig. 5). We have to cover a particle range from 0.08 
m to 0.2m of our dust aggregates with a volume filling 
factor of 0.4. As Wcidcnschilling (1977) assumed com- 
pact particles, we interpolated the curves in his Fig. 5 
with our filling factor of I = 0.4, having in mind that 
the independent variable is the friction time of the dust 
aggregates. Therefore we used his 10-cm curve and ap- 
proximate it with 
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m ( R[AU] Y 
7 \1.1AU ) 



(56) 



To find the particle radius for a given Kepler radius at 
which the drift velocity exceeds 26 m s _1 and therefore 
causes erosion, we combine Eqs. ([53]) and ([56"]) to 



-imp 



m ( R[AU] \ 
~s [lAAU 



0.21m 



set Vi m p to 26 m s 1 



and achieve 



(57) 



(58) 



Eqn. (|58|) is given as a dashed dotted line in the lower 
graphs of Figs. [T0]and[TTJ The particles size above would 
be reduced further in case the impact velocity driven by 
turbulence which is neglected here would be taken into 
account. 

4.1.8. Number Density of Micrometer- Sized Particles 

At Kepler radii larger than 0.2 AU (for a = 10~ 2 ) 
and 0.3 AU (for a — 10~ 4 ), the timescale of the erosion 
is much larger than the coagulation timescale (see top 
graph in Figs. rrUland [rT|) . This means that particles can 
grow to the size predicted in Sect. 14.1.41 and the num- 
ber density of the micrometer-sized particles is given by 
Eq. [32] For Kepler radii smaller than the above values, 
the timescale for the erosion of large particles is much 
shorter than the coagulation timescale. This means that 
particles can only grow as long as their relative motion 
to the gas is smaller than the threshold at which the 
erosion process sets in. To calculate the number den- 
sity of micrometer-sized dust particles produced by the 
erosion process at Kepler radii lower than given above, 
it is assumed that for all particles, which coagulate to 
a size where the relative motion to the gas causes ero- 
sion, the further coagulation is compensated by erosion. 
To get the number density of micrometer-sized particles, 
the production by erosion has to be balanced by their 
re-coagulation. 

In the following, poo denotes the total number den- 
sity at a given Kepler radius if all the dust mass of the 
disk was in micrometer-sized particles and r co is the co- 
agulation timescale to the particle size, v max , at which 
erosion sets in. We assumed for simplicity that this is the 
sa me rate as the coagula t ion to dm-sized particles given 
in iNomura fc N akagawa (200(f) their Fig. 5. However, 
coagulation continues, but all the mass gain of such dust 
aggregates caused by coagulation is immediately eroded 
away in regions where the erosional timescale is faster 



than the coagulation timescale. Thus, we get as the 
source term for monomer dust grains 



dp p _ poo 
dt r rn ' 



(59) 



which gives the production rate of micrometer-sized dust 
particles by erosion. This rate is balanced by the coagu- 
lation of the micrometer-sized particles into dimers, i.e. 



dp p _ o-VBPp 
dt ~ ^2 ' 



(60) 



where Eq. [55] was used to integrate the coagulation, 
which takes place over the whole disk height. Combining 
Eqs. [591 and ITO1 leads to the stationary solution 



Pps — ' 



Poo 



V2 



t co ov b 



(61) 



The number density of the micrometer-sized dust parti- 
cles, resulting from erosion (source term) and coagulation 
(sink term), is shown in Fig. [Tnifor a = 10~ 2 and in Fig. 
[TT] for a = 10~ 4 (second graphs), as a function of dis- 
tance to the central star, respectively. Here, the solid 
line shows the region in which the erosion timescales is 
larger than the coagulation timescale and the dashed- 
dotted line shows where the erosion timescale is smaller 
than the coagulation timescale (see top graphs). We in- 
tegrated all /xm-sized dust particles over the entire disk 
volume and got for their total mass 1.7 x 10 24 kg for 
a = 10~ 2 and 1.1 x 10 24 kg for a = 10~ 4 , respectively. 
Applying a minimum disk mass of 0.0425Mq for the gas 
mass and a total mass ratio between dust and gas of 0.01, 
we get a mass fraction of the ^m-sized dust particles with 
respect to the total dust mass in the disk of 2.0 x 10~ 3 
for a — 10~ 2 and 1.3 x 10~ 3 for a = 10~ 4 , respectively. 

4.2. Consequences for Protoplanetary Disks 

Figs. [10] and [TTJ (top graphs) show that the process 
of impact erosion takes place at orbital distances smaller 
than 5 AU. At orbital distances smaller than 0.2 AU (for 
a = liT 2 ) and 0.3 AU (for a = 10~ 4 ), respectively, the 
erosion timescale of the large dust aggregates becomes 
smaller than the coagulation timescales. Hence, the dust 
aggregates can only grow as long as their relative veloc- 
ity to the gas (an, thus, to the small dust particles) is 
smaller than the threshold velocity for erosion. Dust ag- 
gregates, which grew to larger sizes at larger orbital dis- 
tances, radially drift inward from the outer disk and are 
then eroded to that size. The latter effect is neglected 
in the following, because the number of particles that 
drift in from further out is smaller than the particles 
that are locally produced by coagulation in that area. 
In the second graphs of Figs. [10] and [TTJ the steady- 
state number density of micrometer-sized particles that 
are produced by the above effects and removed by coag- 
ulation are shown. At orbital distances in between ~5 
AU and 0.2 AU (for a = 10~ 2 ) and -5 AU and 0.3 AU 
(for a — 10~ 4 ), respectively, the erosion does not con- 
siderably affect the size of the large agglomerates (but 
still produces many /im-sized dust grains). At orbital 
distances smaller than 0.2 AU (for a = 10 -2 ) and 0.3 
AU (for a = 10~ 4 ), the erosion balances the growth of 
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the large dust aggregates by coagulation. Therefore, the 
number density increase of micrometer-sized dust parti- 
cles towards smaller disk radii is not as steep as further 
out in the disk, because the drift velocity is reduced due 
to the smaller size of the largest particles. Here, radial 
diffusion has been neglected because its timescale is very 
long (see dotted lines in the top graphs of Figs. [TU] and 
[TTj) . In the bottom graphs of Figs. [10] and [11] the con- 
sequences of the erosion on the radii of the largest dust 
aggregates are shown. Note that 70% of the total dust 
mass of the disk is concentr ated in the largest particles 
([Nomura fc Nakagawal 12006). At disk radii larger than 
0.2 AU (for a = 10" 2 ) and 0.3 AU (for a = 10~ 4 ), re- 
spectively, the radii of the particles are nearly unaffected 
by erosion. The growth of dust aggregates at these radii 
is only limited by the inward migration of the particles. 
At smaller disk radii, the growth of the large dust ag- 
gregates is limited by erosion, which means that their 
sizes are reduced to less than 60% (at 0.1 AU) of the size 
without erosion. 

In summary, the process of collisional erosion leads to 
a number density of micrometer-sized dust particles of 
up to 10 7 m~ 3 and limits the maximum dust-aggregate 
size to ~ 0.1 m at orbital distances of 0.3 AU. The ratio 
of /im-sized particles to the total dust mass at a given 
Kepler radius is shown in the third graph of Figs. [TOl 
and [n] for a = 10~ 2 and a = 10 -4 , respectively. For 
a minimum-mass solar nebula model, we integrated over 
the whole disk volume and found a total-mass fraction 
of the /z-sized dust particles of 2.0 x 10~ 3 for a = 10~ 2 
and 1.3 x 10 -3 for a a — 10~ 4 , respectively. 

4.3. Comparison with Observational Data 

To compare the amount of micrometer-sized parti- 
cles fou nd above wit h obse rvational data, we use the 
work of IFurlan et al.l ([2006D . Typical observations of 
iFurlan et al.l ((20061 ) show dust depletion factors of 100 
to 1,000 which means that in their observed disks 10~ 2 
- 10 -3 of the total dust mass still resides in micrometer- 
sized particles, whereas the remainder dust mass is hid- 
den in large dust aggregates. It is noted that the disk 
models used for their data interpretation have a slightly 
shallower density slope with Kepler radius than in the 
minimum-mass disk used here, but we do not think that 
this has an important influence on the observed mass 
depletion factor. In our bimodal model, 2.0xl0 -3 and 
1.3x 10 -3 of the total dust mass is in /Ltm-sized dust par- 
ticles for a — 10~ 2 and a — 10~ 4 , respectively. 70% 
of the total dust mass is contained in large dust aggre- 
gates with a fi lling factor of 40%. This reproduces the 
observations of IFurlan et al.l (|2006l ) very well. 

In case the disks are already evolved and their gas 
density in the midplane at 1 AU has dropped to 10 ~ 9 
kg m -3 , which is 10~ 3 of the value of a disk of mini- 
mum mass with an interstellar dust-to-gas mass ratio, 
the drift velocities of 5 mm-sized dust aggregates (hav- 
ing a filling factor of 40%), corresponding to the drift 
velocities of 2.7 mm com pact particles (s ee above), are 
larger than 26 ms" 1 (see (Wcidcnschilling 1977), Fig. 4) 
and enable erosion much more efficiently. In this case, no 
decimeter- or meter-sized dust aggregates are necessary 
to explain the production and existence of micrometer- 
sized dust grains. This could also explain the cutoff at 



millim eter-sized grains seen in many disks of iRicci et al.l 
((20101) . In our model, the /mi-sized dust particles are pro- 
duced at Kepler radii smaller than ~5 AU. Hence, our 
model predicts that the disk is depleted of micrometer- 
sized dust particles at large Kepler radii while somewhat 
larger grains are still present. 

4.4. Irregular Dust Particles 

As noted above, we used in our experiments monodis- 
perse Si02 spheres as simulated protoplanetary dust par- 
ticles. In reality, protoplanetary dust consists of grains 
with all possible shapes and foll ows a size di stribution. 
As we found in an earlier paper (Blum eTaTl [20061) . the 
tensile strength of dust aggregates consisting of polydis- 
perse irregular Si02 particles is about five times smaller 
than that of the particles we used in the experiment 
described in Sects. [2] and [3] As we can assume that 
the particle-release energy is proportional to the ten- 
sile strength, the erosional threshold should therefore ap- 
proximately occur at 0.5 times smaller velocities. This 
means that the erosional process should work with more 
realistic particles up to disk radii of about 10 A U (using 
the disk model by Schraplcr fc Henninsj (|2004f )). Hence, 
erosion should be properly treated in all simulations in 
which large dust aggregates (exceeding ~10 cm in size) 
play a role and for the entire range of terrestrial-planet 
and the inner part of the giant-planet formation zone. 

5. SUMMARY 

In protoplanetary d isks , realistic dust growth models 
(e.g. IZsom et "al~ll2010f), which do not use a dust ret ention 
mechanism {Okuzumi 2009; IOkuzumi et~aT1l2010allbl) . fail 
to explain the observed large amount of micrometer-sized 
particles (Dullcmo nd fc Dominikl 120051) . Because high 
relative velocities between meter-sized dust aggregates 
and micrometer-sized dust particles exist in protoplane- 
tary disks, we experimentally investigated in a first step 
whether collisions at these relative velocities lead to an 
erosion of the meter-sized aggregates, which could, in 
turn, be a source of micrometer-sized dust particles. We 
found that, after an initial stage in which an impacting 
particle can erode up to 10 dust grains from the fluffy 
dust agglomerate, the agglomerate is surface-compacted 
and passivated against further erosion at velocities lower 
than 26 ms _1 . At larger velocities, the erosion is reduced 
to 10% of its initial strength. In a second step, we ex- 
plained and confirmed the occurrence of impact-induced 
erosion with a numerical model. In a third step, we de- 
veloped an analytical disk model and implemented the 
experimentally found erosive effect. This model shows 
that erosion is a strong source of micrometer-sized dust 
particles in protoplanetary disks. The steady-state solu- 
tion of the model shows that the micrometer-sized dust 
particles encompass a fraction of 0.20% (for a of 10~ 2 ) 
and 0.13% (for 10~ 4 ) of the total dust mass in a disk if 
we assume that 70% of the disk mass is in the form of 
meter-sized dust aggregates with a filling factor of 40%. 
A further outcome of our model is that larger dust aggre- 
gates get eroded very fast at small Kepler radii, because 
they possess higher relative velocities to the micrometer- 
sized particles. As their relative velocities are also de- 
pendent on their size, their growth is limited by ero- 
sion to ~0.1 m at orbital distances smaller than a few 
tenth of an AU. Finally, we used our analytical model to 
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discuss the amount of micrometer-sized particles in the 
observational data of lFurlan et al.l (|2006[ ). We can repro- 
duce their outcome with our erosion model. Assuming 
that the observed disks are already evolved and their 
gas density has dropped to 10 -3 of a disk of minimum 
mass, drift velocities of 5-mm-sized dust aggregates (with 
40% filling factor) are large enough to enable erosion. In 
this case, no decimeter- or meter-sized dust aggregates 



are necessary to explain the production and existence of 
micrometer-sized dust grains. The erosive radius loss of 
grains larger than 5 mm could also explain t he cutoff at 
millim eter-sized grains seen in many disks of iRicci et al.1 

(poioh . 
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